Measuring ultrasmall time delays of light by joint weak measurements 
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We propose to use weak measurements away from the weak-value amplification regime to carry 
out precision measurements of time delays of light. Our scheme is robust to several sources of noise 
that are shown to only limit the relative precision of the measurement. Thus, they do not set a limit 
on the smallest measurable phase shift contrary to standard interferometry and weak-value based 
measurement techniques. Our idea is not restricted to phase-shift measurements and could be used 
to measure other small effects using a similar protocol. 
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Twenty years after the proposal fl)] of Aharonov, Al- 
bert, and Vaidman, several experiments have demon- 
strated the possibility to measure tiny physical effects 
using the so-called weak-value amplification scheme 0- 
0]. These ideas have paved the way for new approaches 
to precision measurements in general, and have triggered 
a great deal of further theoretical and experimental de- 
velopments. 

Recently, Hosten and Kwiat Q were able to experi- 
mentally confirm the spin Hall effect of light by mea- 
suring a polarization-dependent displacement of a laser 
beam to a precision of f A using weak-value amplifi- 
cation. Gorodetski et al. 0] investigated the plasmonic 
spin Hall effect. Dixon et al. Q determined the angle 
of a mirror to a precision of the order of 500 frad by 
measuring deflection of light off the mirror. Several ex- 
periments were proposed in order to enhance the preci- 
sion of the measurement of longitudinal phase shifts of 
light , or amplify the single-photon nonlinearity to a 
measurable effect [6|. An application to charge sensing 
in a solid-state context has also been put forward Q. 
The advantages of weak- value amplification schemes for 
suppressing technical noise were investigated in [1, d, Q , 
and [loj showed how technical noise could even improve 
the precision in this scheme. Ways to optimize the initial 



meter wavefunction were studied in 



fl 



There was a number of attempts to go beyond the 
weak-value formalism. This includes, for instance, 
higher-order expansions for nearly orthogonal pre- and 
post-selected states 12 or the use of full counting 



statistics [17|, and orbital-angular- momentum pointer 
states fisl . The effect of decoherence was investigated in 
Ref. pJl]- The limits of amplification for arbitrary cou- 
pling strength were addressed in 2^ 21 1. Connections 
of the weak-value formalism with the theoretical tools of 
precision metrology were made in [l^, [2^. Weak val- 
ues were also related to quasiprobability distributions of 
incompatible observables [24 1. 

In this Letter, we would like to propose a scheme to 
measure very small time delays, or longitudinal phase 
shifts. The idea is to use weak measurements away from 
the weak- value amplification regime and to exploit the 



full information contained in the correlations induced 
by the time delay between frequency and polarization 
of photons going through the interferometer. This pro- 
cedure is not limited to time-delay measurements and 
could be used for other precision measurements. For ex- 
ample, it is readily applicable to measurements of ultra- 
small beam deflections by slightly modifying the protocol 
of Dixon et al. ^] . The idea of carrying out a full joint 
measurement of two weakly entangled degrees of freedom 
could be relevant in many domains such as charge sens- 
ing in solid state physics precision metrology, and 
gravitational wave detection. 

Following Ref. HI we are interested in measuring with 
high precision a small interaction parameter r which cou- 
ples two physical systems. We will call them system and 
meter as is customary in weak-value protocols, and the 
effect of the coupling is described by an unitary operator 



U 



At Ax 



(1) 



where A is an operator acting on the system and x acts on 
the meter. There are many different situations where the 
measurement of a small coupling constant r is potentially 
relevant. For example the interaction can be an interest- 
ing but very small physical effect, like the spin Hall effect 
of light Q; or T can be a quantity of direct interest, de- 
scribing the angle of a mirror, see [J], or describing the 
amplitude of a gravitational wave in an interferometric 
setup. In this Letter we will consider a setup where a 
time delay in an optical interferometer plays the role of 
the interaction parameter t, as in fsl]. 

The interaction ([1]) with small r is formally similar 
to a weak measurement of the observable A of the sys- 
tem. Therefore by adding a pre- and post-selection of 
the system in states \i) and |/) respectively, see [l[, 
we obtain the weak value of the observable A, A^, = 
(/I A \i) I {f\i) in the meter averages {p) = rReAy,, and 
(x) = — 2r (-01 IV') ImAu,, where jV') is the initial state 
of the meter. Those results are only valid up to lin- 
ear order in r, which is enough for our purposes, but 
the meter averages can also be expressed exactly using a 
generalized form of weak values [25|. Weak- value ampli- 
flcation is based on the fact that it is possible to obtain 
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arbitrarily large weak values by choosing almost orthog- 
onal initial and final states, {f\i) — > 0, thereby amplify- 
ing the small interaction parameter r. The price to pay 
for obtaining large weak values is the small probability 
^ps — of a successful post-selection. Weak- value 

amplification does not increase the statistical information 
because the amplification of the signal is exactly counter- 
balanced by the small success rate of the post-selection. 
However the amplification allows to overcome technical 
limitations, e.g. detector resolution or readout noise. On 
the other hand, stringent post-selections lead to an in- 
creased sensitivity against certain types of noise. For in- 
stance background noise, e.g. stray photons hitting the 
detector in optics experiments, is a limiting factor of how 
small the post-selection can be since the rate of useful, 
post-selected, photons must be significantly larger than 
the background. Another critical type of noise is due to 
errors in the post-selection itself. 

Conceptually, as is the case for all quantum mechanical 
measurements |26j . the crucial information in imaginary 
weak-value amplification is not the average shift of the 
meter, but the correlations induced by the interaction 
between the system and the meter. These are generally 
probed by joint measurements on the meter and system. 
Weak- value amplification constitutes a partial joint weak 
measurement, since measurement outcomes are retained 
only if the system is found in a certain state. Remarkably, 
it turns out that all the useful correlations for the estima- 
tion of T are contained in the small post-selection regime. 
However, additional information on other aspects of the 
experiment is included in the other outcomes. In the 
presence of noise, this additional information turns out to 
be useful to precisely evaluate r. In the following, we will 
present an example where carrying out a full joint weak 
measurement outperforms weak-value amplification. 

Time- delay measurements.- We consider a Mach- 
Zehnder interferometer for laser light, see Fig. [T^. En- 
coding the which-path information in a two-level system 
allows us to write the effect of the time delay as 



yH/2 



(2) 



where ay acts on the which-path space, and H is the 
Hamiltonian of the laser light. The time delay now takes 
the role of the interaction parameter in Eq. ([T]), and the 
photon frequencies (which are contained in H) the role 
of the pointer variable x. The incoming light is evenly 
split into the two arms with opposite circular polariza- 
tions by a polarizing beam splitter. The two arms then 
recombine at another polarizing beam splitter with lin- 
early polarized outputs. The direction of the linear po- 
larization is described by the azimuthal angle (j) on the 
Poincare sphere. The two output ports q = ±1 are mon- 
itored by spectrometers. In the weak-value language, the 
pre-selected polarization state is \i) = (|— ) + |+))/\/2 
and the two post-selected states at detector q are \fq) ~ 



a) PBS (circular) (Ty = -1 PBS (linear, ) 



linear polarizer 



b) 






FIG. 1: (a) Mach-Zehnder interferometer for laser light with 
a time-delay asymmetry r. The incoming beam is linearly 
polarized and then evenly split by a polarizing beam splitter 
(PBS) into the two arms with opposite circular polarizations, 
ay = ±1. The beams are recombined at a second PBS with 
linearly polarized outputs. The direction of the linear polar- 
ization is described by the azimuthal angle (p on the Poincare 
sphere. Two spectrometers labeled by g = ±1 measure the 
spectrum of the outgoing light, (b) Graphical representation 
of the estimation process. The values of r and (j> ^^'^ ex- 
tracted from the spectra measured by the two detectors by a 
maximum-likelihood (ML) estimation technique as described 
in the text. 



+96*'* |+))/\/2. Wc are interested in ultrasmall time 
delays with Wmax''' ^ 1, where w,nax is the maximal fre- 
quency contained in the laser light. 

A laser pulse with normalized spectrum pq{lu), which 
can be interpreted as a probability density, is sent 
through the interferometer. The probability density of 
outcomes is then given by 

Pq{uj;T,(f>) = ipo(w)[l + (7exp(-eV2)cos((/.-wr)] , (3) 

where we have introduced Gaussian fiuctuations of am- 
plitude e of the alignment parameter 0, i.e., convoluted 
the bare probability density with a fluctuation kernel 



1 



27re 



exp 



2e2 



(4) 



that describes fiuctuations around the average alignment 
4>. Note that only the fluctuations with a correlation time 
smaller than the measurement duration are included in e. 
Fluctuations with a longer correlation time will modify 
the effective value of (f). From now on we make the realis- 
tic assumption e ^ 1. It can be shown that fluctuations 
of the alignment of the circular polarizing beam splitter 
and of the linear polarizer can also be encompassed in e. 

Our goal is to find an estimate of the values of the 
parameters t, ip from outcomes of an experiment follow- 
ing the probability distribution ([3]). Although (j) is in 
principle controlled experimentally, its estimation per- 
mits to remove possible systematic errors. The quantity 
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produced by the experiment is a set of observed frequen- 
cies fqiuj), which, in principle, converge to Pq(uj;T,4>) as 
the number N of measured photons is increased. 

We use the maximum likelihood procedure [27| to pro- 
vide unbiased estimates of the values of the parameters, 
see Fig. [T]d. This is achieved by maximizing the log- 
likelihood function l{T,(j)) defined as 



q=±l 



(5) 



Maximizing the log-likelihood yields two equations which 
have to be solved numerically in the general case. How- 
ever, for the special case of almost equal intensities at the 
two detectors, 10 — ujt — §1 ^ 1, analytical expressions 
can be derived 



^|„exp(6V2)^gP,, 



(6) 



^'^eMeV2)(±-j:,Pq{u.)q-j:,Pq) , (7) 



where Pq is the integrated fraction of outcomes in de- 
tector q, and {■)q denotes the average value in detector 
q. The frequency spread Aa; of the initial distribution 

is given by [Aujf = EqP.i^^), " (E, ^(w)?)'- The 
estimates (|6I7[) depend on measurement results and on 
one unknown amplitude e that characterizes alignment 
fluctuations. Remarkably, exp(e^/2) only appears as an 
overall multiplicative factor. Thus, the ultimate error 
Atuh on the estimation of r by assuming e = (since its 
value is unknown, yet realistically e 1) scales with r, 
i.e. only a relative error occurs that does not limit the 
smallest r that can be measured. Equation ([7]) is one of 
the main results of our Letter. 

Aside from the errors due to technical noise, statisti- 
cal uncertainties also contribute to the estimation error. 
For a finite number N of detected photons, the statisti- 
cal errors are provided by the Cramer- Rao bound [27| : 

Arstat > ^/iI-^)rT/N, and A0stat > V(I-i)00/7V, 

where I is the Fisher information matrix given by 



q=±l 



dx pq{x; T, (f)) {dy \0gPq) {dz \0gPq) ■ (8) 



As a side remark, note that the logarithmic derivative 
d-r logPq is connected to the weak value of ay in detector 
q, see Ref. (2^ for a discussion. Asymptotically, i.e., for 
large N, the Cramer-Rao bounds are saturated by the 
maximum-likelihood procedure. In the case of almost 
equal intensities at the two detectors we obtain 



Arstat > 



exp(eV2) 



(9) 



which shows that the fluctuations of (j) do not increase 
significantly statistical noise. The number of detected 



photons N required to obtain a good estimate of r is 
of the order 10/(Allit)^. Estimating a time delay of the 
order of 1 attosecond with ultrashort laser pulses with 
Aw « lO^^Hz would require detecting 10^ photons. A 
typical pulse contains typically 10^"^ photons which would 
be enough to measure time delays of the order of 1 zep- 
tosecond corresponding to a displacement of lOOfm. 

Split detectors.- We would now like to apply this result 
to an actual detector with a finite resolution. This will in 
particular shed light on the roles played by the resolution 
and readout noise of the detector. In principle, the results 
of the previous section could require measuring the full 
distribution of uj without any readout noise. To show 
that this is not the case we consider "split" detectors 
which can only discriminate two spectral regions w > wq, 
leading to a measurement result r = and u < ujq, 
leading to r — —1. We also add readout noise from the 
outset. We will see that our conclusion from Eq. ([T]) (viz. 
that there is no absolute lower limit on the smallest value 
of T that can be measured) survives even in this extreme 
case. 

To make up for losing the detector resolution some 
a priori knowledge of the initial frequency distribution 
Pq{(jj) is required. For analytical calculations we will as- 
sume an initial Gaussian distribution, but the results 
will depend only quantitatively on the shape. In prac- 
tice the distribution should be measured and the cal- 
culations done numerically. We thus assume po{uj) w 
g-(i^-i^o)V2(Ai^)^^^/2^^^ where the tails should, in prin- 
ciple, be truncated since < oj < Wmax- The probabilities 
of the two possible outcomes r = ±1 at the two detectors 
q = ±1 at second order in the small quantity Awr read 

Prq = i[l + <7exp(-eV2)(l- ^(Awr)2)cos(/)] 



+rqexp{-e'^/2) 



AUJT 



2^27r(l + (n/Auj)^) 



sin0, (10) 



where we allowed for a frequency detection uncertainty 
of order modeled as Gaussian white noise. 

Denoting the measured probabilities by frq, the max- 
imum likelihood estimation can be analytically carried 
out in the regime |sin(/)| 3> Alot, and yields coscj) — 
exp(e^/2)(P+ — P_) in terms of the integrated fraction 
of outcomes Pq = frq, and 



y/2n{l + jn/Aujy) 
8AL^v/exp(-e2) _ (p+ _ p_)2 ^ 



rq.f, 



rq 



P-) 
(11) 

The first terms in the expansion of r in the fluctuations 
e and detection noise VL can now be calculated, and we 
finally obtain 



T = To 



1 



sm( 



1 / VL 

2 V A^ 



(12) 
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where tq is the estimated value of t in the absence of 
noise, see Eq. (fTTj) with e and ft set to zero. Again the 
noise does not set an absolute limit on the precision of 
the estimation of r but only a relative precision. The 
frequency spread Aw reduces the effect of readout noise. 
Moreover we observe that in order to minimize the ef- 
fect of fluctuations of </), we have to work in the regime 
sin (/) « 1 and not in the weak- value amplification regime 
of simp « where the effect of fluctuations increased. 

The statistical uncertainty is given by the Fisher in- 
formation matrix, which is diagonal in this case, through 
the Cramer-Rao bound 



At,, 



> 



27r 



4AwViV 



1 + ::: 



(13) 

Hence the considered fluctuations do not jeopardize the 
estimation of r away from the weak-value ampliflcation 
regime. We also note that using split detectors leads to a 
modest increase of statistical noise by a factor •\/27r with 
respect to Eq. ©. Equations (|12ll3p demonstrate that 
even for low-resolution detectors our scheme is robust 
against readout noise and alignment errors. 

Comparison to existing schemes - Standard interfer- 
ometry compares two probabilities given by the sine and 
cosine of the total phase shift given by the combination 
(j) — LOT. This leads to two difficulties: first, to estimate 
r precisely, the laser frequency w has to be highly sta- 
bilized. Secondly, the alignment <j) cannot be separated 
from the effect of r, i.e., alignment errors are a limiting 
factor to the precision achievable having complete statis- 
tical information. This ultimate precision, which cannot 
be increased by acquiring more measured data, is given 
by Atuu = e/uj. In the procedure proposed by Brun- 
ner and Simon 0j that uses the imaginary part of the 
weak value, the first issue is solved since a large frequency 
spread Acj is advantageous for the precise evaluation of 
T, which is also true in our scheme. However, the sec- 
ond issue is only partially addressed in alignments 
errors are still a limiting factor, Atuh = Ce, where the 
proportionality constant C ~ 10~^* s is three orders of 
magnitude smaller than for standard interferometry. 

In contrast to that, a major advantage of our scheme is 
to remove systematic errors as well as fluctuations of the 
alignment as a limiting factor to the ultimate precision 
of the time-delay measurement. Alignment fluctuations 
lead to a relative error Aruit = e^r/2 on the estimation of 
r, see the discussion after Eq. (O and the illustration in 
Fig. [2l This is made possible by working away from the 
weak-value amplification regime and using all the infor- 
mation contained in the correlations between frequency 
and polarization of the photons to perform a simultane- 
ous estimation of 4> and r. 

Finally, we would like to mention that if the detector 
can only measure events at a small rate due to detector 
saturation, small post-selection probabilities as realized 




t[10"^^s] 



FIG. 2: Ultimate precision limit Atuh on the estimation of 
the physical time delay r as a function of r, in the presence 
of fluctuations of the alignment (f), see Eq. Q. We assume 
e = 0.02 rad for the plot. The dashed line shows the constant 
(zeroth-order) contribution of fluctuations to Atuu = Ce for 
weak-value amplification, here we assume C ~ 0.25 x 10~^* s 
like in The solid line is the result of our scheme, Atuu = 
e^r/2. For r < 2C/e (shaded area), our scheme outperforms 
weak-value amplification . 



in weak- value amplification schemes permit to effectively 
increase the rate of measurements [8| . This allows to ob- 
tain better statistics, however, if the measurement is lim- 
ited not by statistics but for other reasons, our method 
appears preferable. 

In conclusion, we have proposed a technique that could 
be useful to determine very small time delays, or longi- 
tudinal phase shifts, due to its robustness with respect 
to various noise sources. The time delay induces correla- 
tions between frequency and polarization of photons go- 
ing through the interferometer, and our scheme exploits 
the full information contained in these correlations. The 
idea is not limited to time-delay measurements and could 
be used for other precision measurements. In a broader 
setting, the idea of carrying out a full joint measurement 
of two weakly entangled degrees of freedom could find 
applications in many domains such as charge sensing in 
solid state physics, precision metrology, and gravitational 
wave detection. 
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